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ABSTRACT 


Context. The 340-year old supernova remnant Cassiopeia A, located at 3.4 kpc distance, is the best-studied young core-collapse 
supernova remnant. Nucleosynthesis yields in radioactive isotopes have been studied with different methods, in particular, methods 
for production and ejection of 44 Ti and 56 Ni, which originate in the innermost regions of the supernova. 44 Ti was first discovered in 
this remnant, but is not seen consistently in other core-collapse sources. 

Aims. We aim to measure radioactive ^Ti ejected in Cassiopeia A and to place constraints on velocities of these ejecta by determining 
X- and y -ray line-shape parameters of the emission lines. 

Methods. We analyzed the observations made with the SPI spectrometer on INTEGRAL together with an improved instrumental 
background method, to achieve a high spectroscopic resolution that enables interpretation for a velocity constraint on 44 Ti ejecta from 
the 1.157 MeV y-ray line of the 44 Sc decay. 

Results. We observe both the hard X-ray line at 78 keV and the y-ray line at 1157 keV from the 44 Ti decay chain at a combined 
significance of 3.8 cr. Measured fluxes are (2.1 ± 0.4) 1(L 5 ph citT 2 s -1 and (3.5 ± 1.2) 1(L 5 ph cnT 2 s _1 , which corresponds to 
(1.5 ± 0.4) 1(L 4 and (2.4 ± 0.9) 1(L 4 M 0 of 44 Ti, respectively. The measured Doppler broadening of the lines implies expansion 
velocities of 4300 and 2200 km s _1 , respectively. By combining our results with previous studies, we determine a more precise 
estimate of ejected 44 Ti of (1.37 ± 0.19) 1(L 4 M 0 . 

Conclusions. The measurements of the two lines are consistent with previous studies. The flux in the line originating from excited 
44 Ca is significantly higher than the flux determined in the lines from 44 Sc. Cosmic-ray acceleration within the supernova remnant 
may be responsible for an additional contribution to this line from nuclear de-excitation following energetic particle collisions in the 
remnant and swept-up material. 

Key words. (Stars:) supernovae: individual: Cas A - Nuclear reactions, nucleosynthesis, abundances - Stars: massive - ISM: super¬ 
nova remnants - Gamma rays: ISM - Techniques: spectroscopic 


1. Introduction 

Cassiopeia A (Cas A) is the closest young remnant of a core¬ 
collapse supernova and therefore a prominent study object in 
many wavelength regimes. It may have been seen in AD 1680 
by Flamsted (Ashworth 1980 Green|2002]), b ut was not as brigh t 


as expected from its proximity of 3.4 ': ? kpc (Reed et al. 


1995). 


Apparently, the supernova itself was occulted by interstellar dust 
and gas, and its brightness and 56 Ni production are therefore only 
inferred indirectly (Erik sen et al.|2009]). The date of the event is 
not precisely settled; |Thorstensen et al.j ( j200lj ) inferred the su¬ 
pernova event at AD 1671.3+0.9 by measuring the proper mo¬ 
tion of the supernova remnant (SNR) bright main shell as well 
as faint, high-velocity, outer ejecta knots, and assuming free ex¬ 
pansion. 

Cas A was the first supernova in which the characteristic 
line at 1157 keV from 44 Ti decay had been reported (Iyudin 
et al. 1994); this line arises from de-excitation of 44 Ca, the 
final daughter nucleus in the decay chain 44 Ti—> 44 Sc—> 44 Ca. 
Later measurements focused on the equally bright lines from 
the first stage in the 44 Ti decay chain, the de-excitation of 44 Sc. 
In particular, the studies with INTEGRAL/IBIS (Renaud et al. 
|2006j) and NuStar (Grefenstette et akj|2014| settled the amount 
oP^Ti to — 1.4 10 4 M 0 (IBIS: lV“f 1(T 4 M 0 ; NuStar: 
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(1.25+0.3) 1CL 4 M 0 ). This is significantly higher than the yields 
predicted for typical supernova models, which suggested that 
Cas A represents a rare subclass of 44 Ti ejecting supernovae, and 
expl osion asymmetries have been thought t o cause such anoma¬ 


lies (Nagataki et al. 1998, The et al. 2006 


The 44 Ti image re¬ 
cently obtained by the NuStar hard X-ray telescope ( |Grefenstette| 


ly obtair 

et al. 2014) has shown that Cas A ejecta carrying 44 Ti appear 
in clumps, re-affirming that Cas A did not explode as a spheri¬ 
cally symmetric supernova. Other observations had shown this 
before from X-ray measurements, for example, and especially 
from optical/IR light echo spectra, revealing this asymmetry in 
the supernova photosphere. Furthermore, it was shown by these 


light echoes that Cas A exploded as a type lib supernova (Krause 
let al.||2005 |20081 |Rest et al.||2008| |201 1| see also discussion 
by Wheeler et al. ( |2008| )). 


Earlier gamma-ray spectroscopy with INTEGRAL/SPI had 
provided some velocity constraints on the radioactive ejecta 
from the inner supernova, with a lower limit on 44 Ti ejecta of 
500 km s 1 from a nondetection of the high-energy line at 1.157 
MeV, which was attributed to Doppler broadening (Martin et al. 


2007). Doppler broadening is proportional to photon energy; 
therefore, a background dominated measurement such as with 
SPI would be more sensitive to a low-energy line, while the 
high-energy line would drown in the background, especially if 
part of the 44 Ti ejecta were moving at high velocities. NuStar 
detectors also provide a sufficiently high spectral resolution for 
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velocity information. Grefenstette et al. ( j2014| ) reported an over¬ 
all velocity of fastest 44 Ti clumps of (5350 ± 1610) km s 1 , 
simultaneous to a bulk line-of-sight velocity in the range of 
2000 km s _1 . The NuStar result is only based on the detection of 
the 68 keV line, although both hard X-ray lines have been mea¬ 
sured. Our recently improved method for handling instrumen¬ 
tal background suggested re-analysis of INTEGRAL data from 
Cas A with the aim to revisit and refine velocity constraints from 
high-resolution spectroscopy with the SPI instrument of all three 
lines resulting from the 44 Ti decay chain. 


44 Ti and 56 Ni arise from nuclear burning in the innermost 
parts of supernovae very close to the mass-cut separating ejecta 


from the compact remnant (Diehl & Timmes |1998 1 . Therefore 
44 Ti and 56 Ni are expected to also be co-spatial when they de¬ 
cay i n the expanding SNR (|The et al. |1998| Hoffman et al. 
2010[ Magkotsios et a l.[|2010| ).~But as this inner region is prob¬ 
ably characterized by dynamical instabilities and simultaneous 
inflows and outflows of material, it remains unclear how and 
where this separation between material accreting onto the com¬ 
pact remnant and the ejecta occurs, and thus how much of the 
44 Ti can end up in the ejecta (Fryer et al. 2008 Wongwathanarat 


et al. 2013} Popov et al.|2014j>.' 



Galactic Longitude 


Fig. 1 . Exposure map of the sky with the SPI telescope on 
INTEGRAL for the data used in this analysis (Mar 2003 - 
Jan 2014). Cas A is marked with a star at its position (l/b) = 
(111.735°/ - 2.130°). The total exposure time in the PCFOV is 
10.84 Ms, while Cas A was seen fully coded within an exposure 
of 5.15 Ms. 


2. Data and analysis 


56 Ni radioactive decay occurs with a first decay to 56 Co after 
- 8 days (da Cruz et al']|1992 1 , that is, when supernovae are 


2.1. Instrument, exposure, and general analysis method 


still sufficiently dense to absorb even gamma-rays at MeV ener¬ 
gies and convert this radioactivity energy into thermal emission 
which makes supernovae shine at UV to IR wavelengths (|Isern 


eTar|[2008] [Ropke et al.|20121 >. As the supernova envelope be¬ 


comes more transparent, fewer gamma-rays thermalize, in par¬ 
ticular when the second decay stage from 56 Co to 60 Fe at r ^ 111 
days ( |Funck et al.|1992 1 produces 0.847 and 1.238 MeV gamma- 


the explosion (Isern et al. 2008). 44 Ti has a considerably longer 
decay lifetime, and thus gamma-ra vs will escape readi ly. 44 Ti 
decays to 44 


et al. 

1998 

Norman et al. 1998 1 Wietfeldt et al. 1999; Hashimoto 

et al. 

2001 

Ahmad et al. 2006). In this first decay stage, gamma- 


rays of 67.87 and 78.36 keV are emitted. The subsequent decay 
of 44 Sc to 44 Ca occurs after r - 5.73 hours (Audi et al. 2003 


producing a characteristic gamma-ray line at 1157.02 keV. But 
in spite of this long decay time, significant deposit of radioac¬ 
tivity energy is expected in the late SNR from positrons emitted 
during these /3 + decays. 

For SN1987A, this has been investigated: The study of the 
late light-curve i n SN1987A (|Lundqvist et ah] 2001 Fransson 
|& Kozma||2002[ |Jerkstrand et al.||2011[ ) obtained an inferred 
44 Ti mass from optical and infrared spectra of SN1987A to 
(1.5+0.5) 10 4 M 0 . The 44 Ti that mainly powers the late light- 
curve (years after the explosion) has a direct effect on the to¬ 
tal absolute flux level of the optical-to-infrared spectra, but also 
depends on the estimated extinction. Corresponding gamma-ray 
emission has recently been detected by Grebenev et al. (2012 1 
with IBIS; this direct measurement constrained the 44 Ti mass 
of SN1987A to (3.1+0.8) 10~ 4 M 0 (see, however, Seitenzahl 
et al. (2014), who obtained (0.55+0.17) 10 -4 M Q from a de¬ 


tailed, multi-component light-curve model). Beyond SN1987A, 
the Cas A SNR is the only other core-collapse event where this 
consistency between nucleosynthesis radioactive elements and 
observed emission can be studied. 


rays. Therefore, these characteristic gamma-rays are expected to 
be observable with gamma -ray telescopes a few months after 


The INTEGRAL space gamma-ray observatory (Winkler et al. 
2003) carries the spectrometer instrument SPI as one of its two 
main instruments, measuring gamma-rays in the energy range 
of 20 keV to 8 MeV with a spectral resolution of ~ 2.3 keV at 
1 MeV ( Vedrenne et al.|2003j|Roques et al.|2003| ). SPI features a 
camera consisting of 19 Ge detectors with high spectral resolu¬ 
tion, which measures celestial gamma-rays through shadograms 
imprinted by a tungsten coded mask 171 cm above the Ge cam¬ 
era. A high instrumental background from cosmic-ray activation 
of spacecraft and instrument materials makes it challenging to 
distinguish these shadowgrams from the total signal. SPI data 
consist of energy-binned spectra for each of the 19 Ge detec¬ 
tors of the telescope camera, typically taken in 30 min exposures 
dithered in ~2° steps in a certain sky region around the target 
of interest. For our analysis, we used exposures in which Cas 
A was at least in the partially coded field of view (PCFOV) of 
SPI of 34° x 34° accumulated over almost eleven years of the 
INTEGRAL mission. Solar flares mainly affect the background 
rate and are difficult to model. We therefore excluded data of any 
orbit in which a solar flare occurred until the overall background 
returned to normal. Perigee passages around Earth, involving 
a higher radiation dose when the satellite is passing the Van 
Allen radiation belts, additionally affect the background rate, 
and we consequently chose data from orbital phases [0.15,0.85]. 
Background variations due to cooling (heating) after (before) an¬ 
nealing periods were removed accordingly. Finally, our 11-year 
observation set consists of 5604 instrument pointings, resulting 
in a total exposure time of 10.84 Ms (PCFOV) in which Cas 
A was seen fully coded by the SPI coded mask aperture during 
5.15 Ms (see Fig. |T](. 

In our spectroscopic analysis we generally fit scaling fac¬ 
tors of an adopted model of the sky intensity distribution plus 
a model of instrumental background to the set of spectra accu¬ 
mulated during the multiyear observations from our 19 Ge de¬ 
tectors. This resembles linear regression of specific, given time 
sequences, incorporating the instrumental (imaging) response 
function (IRF) for the sky contribution to the data. Multiple com¬ 
ponents of sky and background can be used as required by ade- 
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quacy of the fit (mainly background) and the scientific question 
(sky components). Hence, data d k of energy bin k are modeled 
as a linear combination of the Ni sky model components M l j as 
convolved with the instrument response matrix Rj k , and the Nb 
background components Bj k : 

_ N, N,+N b 

dk = K j k ^ OiMij + ^ 8jBj k . ( 1 ) 

j i=l i=JVi+l 

The comparison between scaled models and measurements is 
performed in data space, which consists of the counts per energy 
bin k measured in each of the SPI detectors j for each single tele¬ 
scope pointing as part of the complete observation. Our adopted 
sky intensity distribution for Cas A is a point source, given the 
SPI imaging resolution of 2.7 degrees and the total SNR diam¬ 
eter of ~5 arcmin of Cas A. The model for the behavior of the 
instrumental background is derived from a separate analysis us¬ 
ing all available data during and around the time of Cas A expo¬ 
sure. We repeat this fit (model scaling) for 0.5 keV wide energy 
intervals to obtain the source intensity spectrum, that is, without 
assuming a prior spectral shape. 

A priori, we do not assume the processes that lead to the 
low-energy lines (68 keV and 78 keV) that originate in excited 
44 Sc, to be the same as for the high-energy line (1157 keV) 
from excited 44 Ca. For example, cosmic-ray (CR) collisions with 
ambient 44 Ca at several tens of MeV/nucleon effectively pro¬ 
duce nuclear excitation, followed by de-excitation photons at 
1.157 MeV. Because 44 Sc is short lived, no corresponding excita¬ 
tion of 44 Sc will occur through this cosmic-ray process, meaning 
that no additional emission is produced in the 68 and/or 78 keV 
lines. We therefore treat the hard X-ray lines and the gamma-ray 
line indepedently. 


2.2. Background modeling 


In earlier SPI analyses, the detailed background line informa¬ 
tion was neglected or used only indirectly as background trac¬ 
ers, obtained from cosmic-ray intensity variations or from adja¬ 
cent energy bands, for example. These were separately adjusted 
by a set of fit parameters for each energy bin in the measured 
spectra (Wang et al. 2007 2009} . To improve the spectroscopic 
sensitivity compared to previous SPI analysis results, we im¬ 
plemented a new approach for modeling the instrumental back¬ 
ground. 

We now investigate the spectroscopic signatures of nuclear 
reaction physics that occur in the instrument and satellite in great 
detail by separately tracking the shape and intensity variations 
of characteristic lines in each Ge detector. First, a high-precision 
cumulative spectrum is used in the two energy bands, ranging 
from 60 keV to 85 keV, and from 1143 keV to 1175 keV, re¬ 
spectively, to also identify the weak instrumental line signatures. 
Then, a three-day integration is used to obtain sufficient statis¬ 
tical precision to determine spectral shape parameters of each 
identified background component for all detectors by fitting the 
spectra using a Metropolis-Hastings algorithm in a Monte Carlo 
Markov chain. We then fix the shapes of the spectral templates 
during one orbit of the selected data set and scale the ampli¬ 
tude of each component for the entire Ge camera by pointing 
to properly trace intensity variations of each background com¬ 
ponent. Long-term investigations showed that the relative inten¬ 
sities of a particular spectral background feature present in one 
detector with respect to the mean intensity of the same feature 
in all detectors (detector ratios) are constant in time and thus are 
maximally independent of possible celestial signals ( [Diehl et ah] 


2014 


2015 ” Therefore, the spectral shape and detector ratio 


constancy that is expected for each individual background com¬ 
ponent is imprinted onto the short-term (pointing-to-pointing) 
variations, while the mean amplitude (over all detectors) of a 
particular spectral feature is still allowed to vary to normalize 
the predicted background count rate. 

In principle, each background component may have its own 
properties, and in particular, detector-to-detector intensity ratio, 
but we clearly identify physically plausible classes: Ge back¬ 
ground lines provide a higher count rate in the inner Ge detec¬ 
tor elements which are surrounded by Ge, while Bi background 
lines provide a lower count rate in the inner Ge detector elements 
because the BGO anticoincidence shield surrounds the Ge cam¬ 
era. For details see Siegert et al., in preparation (2015). 

The instrument spectral response is represented as a pa¬ 
rameterized semi-analytic function, that is, the convolution of 
a Gaussian with an exponential tail toward lower energies to de¬ 
scribe the line shape. The spectral model we use in our line fit¬ 
ting consists of a linear continuum 


Cj(E\ Cqj, Cij) - Coj + C i j ■ E 


( 2 ) 


and a set of instrumental lines at positions Eqjj , which are 
needed to characterize the spectrum. Each Gaussian line i 


Gij(E, Aojj, Eqjj, cr jj) — Aojj ■ exp 


(E - E 0 , u ) 

2-cr? 


2 \ 


(3) 


is convolved with an exponential tail function 


Tij(E\ t ij) = — exp 



(4) 


that describes the impact of cosmic radiation that gradually de¬ 
teriorates the charge collection efficiency of each detector j. The 
raw spectra of SPI are background dominated; the instrumental 
background is more than 99% of the total, with a small contri¬ 
bution from a celestial signal. Fitting the raw spectra for each 
detector and a three-day cumulative sample results in a very pre¬ 
cise spectrum of what is observed from instrumental background 
gamma-ray lines and continuum. This ’’continuum” includes the 
instrumental background continuum (> 99%) as well as any pos¬ 
sible continuum emission from celestial sources. We are mainly 
interested in the emission of gamma-ray lines; therefore, this 
approach largely also suppresses the continuum from celestial 
sources, and any resulting offset from zero in the finally derived 
spectra can be considered as a measure of systematics of this 
continuum treatment. For each energy bin, a discrete background 
detector pattern (in time and position in the sky) is predicted, 
which serves as the only background component Bj k in Eq. O- 
The degradation of Ge detectors from cosmic-ray irradiation and 
their restoration in annealings results in a time-variable width 
and asymmetry of the spectral response. This variation domi¬ 
nates all other spectral changes and is found consistent across 
the SPI energy range. 


3. Results and discussion 

3.1. 44 Ti gamma-ray line measurements 

Our analysis shows the expected signature of a (Gaussian- 
shaped) line at the expected energies for the two independent 

1 This is true for the time during one camera configuration. Within 
these 11 years of data, 4 detectors out of 19 have failed in flight, thus 
leading to a new camera configuration - in total, five different configu¬ 
rations. 
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spectra. The significances of the lines are 2.2 cr (1157 keV) 
and 3.1 cr (78 keV), with a combined 44 Ti detection signifi¬ 
cance of 3.8 cr. Unfortunately, systematic effects due to strong 
background features below 70 keV cannot be reliably modeled; 
therefore no constraint on the 68 keV line could be obtained. 
The mean reduced x 2 of the derived fit per energy bin is 1.06 
( x 2 = 98241, d.o.f. = 92517), suggesting an additional system¬ 
atic uncertainty of - 6 %. This is about eight times lower than the 
statistical uncertainties. 

In Figs. [2] and [3] we show the derived spectra around the 
1157 keV line and the 78 keV line. The spectrum in Fig. [2] 
is well represented by a constant plus a Gaussian that is cen¬ 
tered at (1158 ± 3.6) keV and broadened with respect to the 
instrumental resolution. We obtain a total flux of F 1157 = 



Energy [keV] 

Fig. 2. Spectrum obtained from a source at the position of Cas A 
at 1157 keV (black crosses, 2.5 keV energy bins). The spectrum 
has been fitted by a constant plus a Gaussian that is centered at 
(1158 ± 3.6) keV and broadened (8.94 ± 1.42) keV (FWHM). 
The derived expansion velocity of the 44 Ti ejecta is (2200 ± 
400) km s' 1 . The measured flux is (3.5 ± 1.2) 10' 5 ph cm' 2 
s' 1 , which can be converted to an observable 44 Ti mass of 
(2.4 ± 0.9) 10' 4 Mg at the time of the explosion. The cumula¬ 
tive long-time spectrum (mainly background) is shown as a gray 
histogram. The laboratory-determined energy of the excited 44 Ca 
line is marked with an arrow. 

(3.5 ± 1.2) 10~ 5 ph cm' 2 s' 1 from this fit. This can be trans¬ 
lated into an observable 44 Ti mass at the time of the explosion 
(t = 0) of (2.4 ± 0.9) 10' 4 Mo, by 

Mu Ti (t = 0) = And 1 F l 157 (f) ■ 44m„ • Tu Ti ■ e\p(t/Ta Ti )/W n 5T, (5) 

taking into account that the characteristic lifetime of 44 Sc 
(5.73 h) is much shorter than of 44 Ti (86 a). In Eq. ([ 5 ]), d = 
3.4 4 JJ- 2 kpc is the distance to the SNR, t = (337.2 ± 6.4) a its age, 
incorporating the eleven-year observation time span, F\ 157(f) the 
measured flux in the 1157 keV line today, m„ the atomic mass 
unit, T 4477 = (86 ± 0.5) a the characteristic lifetime of the 44 Ti 
decay, and VVj [57 = 99.9% the respective branching ratio. 

The uncertainty on this 44 Ti mass value is mainly due to 
the distance to Cas A and the measured flux. In addition, the 
44 Ti half-life time has an impact on the derived mass; its uncer¬ 
tainty is taken into account above. Because statistical fluctua¬ 
tions in the measured spectra are large, the fit and the resulting 
mass might suffer from systematic effects in the assumed spec¬ 
tral shape (constant plus Gaussian); we estimate this effect by 



70 75 80 85 

Energy [keV] 


Fig. 3. Same as Fig. [ 2 ] but at 78 keV (1.5 keV energy bins). 
The best-fit values are (78.9 ± 1.5) keV for the centroid, 
(1.9+ 0.6) keV for the width (FWHM), yielding an expansion 
velocity of (4300 ± 1600) km s' 1 , and (2.1 ± 0.4) 10 ' 5 ph cm ' 2 
s ' 1 for the flux, which is equivalent to (1.5±0.4) 10 ' 4 M G of 44 Ti. 
The raw count spectrum (gray) increases by two orders of mag¬ 
nitude below 70 keV with respect to 80 keV, which introduces 
strong systematic effects. 


determining an upper limit of the constant offset, which would 
then reduce the derived mass by about 33%. 

The high-energy 44 Ti decay line has a FWHM of (8.94 + 
1.42) keV and is thus broadened with respect to the instrumen¬ 
tal resolution at 1.157 MeV of ^ 2.4 keV ( |Attie et ak1|2003| . 
This FWHM can be converted to an estimate of the expan¬ 
sion velocity of the 44 Ti ejecta of (2200 ± 400) km s' 1 . This 
value is significantly lower than values reported in recent pa¬ 
pers (Renaud et al. 2006 Grefenstette et al. |2014[ >. The line 
centroid is found as (1158.0 ± 3.6) keV and thus agrees with 
the laboratory-determined energy, suggesting only minor bulk 
average motion of 44 Ti in Cas A relative to the solar system, 
(vHuik = (-250 ± 1000) km s' 1 ). 


For the 78 keV line (Fig. [3]), we obtain a line flux of 
F 78 = (2.1 ± 0.4) 10 ' 5 ph cm ' 2 s' 1 . Using Eq. (JsJ and replac¬ 
ing W 1157 with W 78 = 96%, F 78 is converted to a 44 Ti mass of 
(1.5 ± 0.4) 10 ' 4 M g . Within uncertainties, the two derived 44 Ti 
masses are consistent. The 78 keV line is not significantly shifted 
either, confirming absence of average bulk motion. The mea¬ 
sured Doppler broadening of (1.9+ 0.6) keV (FWHM) yields an 
expansion velocity of (4300 ± 1600) km s ' 1 for the 44 Ti ejecta. 
This value is consistent with the result from the 1157 keV line 
and also with the value measured by NuStar. The apparent dif¬ 
ferences between values derived for the 78 and 1157 keV lines 
reflect the energy dependence of the resolving power of the SPI 
instrument: While at 78.4 keV, the resolving power of the spec¬ 
trometer is /?78 = 78.4/1.6 = 49, it is R 1157 = 1157/2.4 = 482 at 
1157 keV, meaning that SPI derives better line width constraints 
at higher energies; this adds to the linear energy dependence of 
the Doppler effect. Figure [4] shows the 2cr upper limit on the 
expansion velocity assuming the best-fit values for flux and cen¬ 
troid for the two detected lines. While the constraint from the 
78 keV line is rather flat with a shallow minimum (gray line 
in the figure), the 1157 keV line (black line) provides a well- 
defined minimum or constraint. 

Representing the line shapes by Gaussians is an empiri¬ 
cal approach. Alternatively, a physical model can be adapted: 
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68 keV line ( 44 Ti(EC) 44 Sc') 

■ 78 keV line ( 44 Ti(EC) 44 Sc) 

■ 1157 keV line ( 44 Sc(EC) 44 Ca') 



0 


OSSE BeppoSAX ISGRI NuStar SPI 
COMPTEL (This work) 


Fig. 4. Ax 1 vs. expansion velocity for the two detected lines, in¬ 
dividually. The values are calculated assuming the best-fit val¬ 
ues in each case, only letting the Gaussian width vary. Because 
of the higher resolving power of SPI at 1 MeV with respect to 
lower energies, the value for the expansion velocity derived form 
the high-energy line (solid black) is more constrained than that 
from the low-energy line (solid gray). For comparison, the 2<x 
limits are indicated by the thick dashed line. 


Assuming a homogeneously expanding thin shell where the shell 
radius is much larger than the distance to the observer, the ex¬ 
pected ideal line shape can be described by a tophat function 


S(E-,F 0 ,Eo,R) = F 0 -Q /? 


1 


2 RE f) 


( 6 ) 


where Fq is the measured flux, Eq the bulk motion, R the ex¬ 
pansion velocity, and c the speed of light (Kretschmer, K. A. 
201 1|>. This line shape has to be convolved with the spectral re¬ 


sponse function of SPI (see Eqs. ([3]) and (|4])). Fitting the spec¬ 
tra using Eq. leads to lower velocities in both cases, that 
is, (3900 ± 800) kms- 1 (78 keV) and (1500 ± 400) krns' 1 
(1157 keV), respectively. The measured expansion velocities, 
however, are similar to those from a Gaussian fit, and system¬ 
atic differences are smaller than uncertainties from statistics. 

When the results of |The et al.| ( |1996) obtained w ith OSSE 
on CGRO with the COMPTEL result by Iyudin et al. ( 1997) , the 


BeppoSAX result (Vink et al. 


suit (Renaud et al. 


2001) , the INTEGR AL/IBIS re- 


2006), and the NuStar result (Grefe nstette| 


et al. 2014l are combined with our work through a weighted 
mean (weighted by the inverse variance of each measurement), 
the 44 Ti mass is (1.37 ± 0.19) 10~ 4 M 0 . 


3.2. Flux consistency check 

We note that the high-energy line appears to yield a systemat¬ 
ically higher 44 Ti mass, both from our analysis and in Iyudin 
et aI7|( 1994) and Iyudin et al. (1997 1 , with respect to the hard X- 
ray line based 44 Ti mass determinations. Our derived 44 Ti mass 
from the 78 keV line is consistent with other published studies 
that measured the low-energy lines. The derived 44 Ti mass per 
instrument is shown in Fig. [5] together with the weighted mean 
of all measurements (hatched area). This systematic difference in 
44 Ti flux values becomes even more obvious when the measure¬ 
ments of the 44 Sc and 44 Ca lines, taken with different instruments 
in the past two decades, are compared and plotted on a time axis 


Fig. 5. Derived 44 Ti masses of different instruments. All data 
bars are drawn from the measurements with only the given in¬ 
strument (no combinations), including the uncertainties in half- 
life time and explosion date. The OSSE and COMPTEL values 
are corrected for the distance to Cas A. Upper limits (if given) 
are shown as dashed lines. The yields measured by COMPTEL 
and SPI using the 1157 keV line are systematically higher than 
values using the low-energy lines. The calculated weighted mean 
of all shown measurements is illustrated as the hatched area and 
is (1.37 + 0.19) 1(T 4 M 0 . 



Year 


Fig. 6. Compilation of y-ray flux measurements in Cas A from 
different instruments during the years 1992 to 2014 (see text). 
To each family of lines (low- and high-energy), an exponential 
decay function has been fitted, fixing the half-life time of 44 Ti 
to 86 a and the explosion date to AD 1671 (solid lines). The lcr 
uncertainties are separately shown by dashed gray lines for each 
fit. The discrepancy can be expressed as a constant flux level of 
(2.29 ± 0.62) 10^ 5 ph cnT 2 s _1 . 


(Fig. |6). The values in Fig. |6| are taken from |The et ah] (|1996 


Table 3) with OSSE, |Iyudin et al.| ( fl997> with CO MPTEL, jVink 


|etal.| ( |200T] Table 1) with BeppoSAX, |Renaud et al.|(|2006[ Table 
1) with IBIS, this work with SPI, and |Grefenstette et al.| ( |2014| 
Table ED2) with NuStar (from left to right). These values only 
consider the value measured with the mentioned instrument it¬ 
self, hence no double-counting of flux values. The 1157 keV flux 
values ( |Iyudi n et al. (1997 1 and our study) are found to be sig¬ 
nificantly higher than the low-energy line flux measurements. 
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We investigated this systematic by treating the measurements 
of the hard X-ray lines independently of the 1157 keV gamma- 
ray line. We performed an F-test to check whether the derived 
spectra can be represented by only one flux value (case I) or if 
the preference is for two values (case II). The corresponding x 2 
values are ;^ 2 = 15.14 (16 dof), and x]i = 12-27 (15 dof), re¬ 
spectively. The F-value of this test is F = 3.51, for which the F- 
statistic gives a probability of ~ 1%, corresponding to 2.5 cr, that 
the two spectra cannot be represented by the same flux value. 
Alternatively, we fixed the explosion date to AD 1671 and the 
characteristic lifetime of 44 Ti to 86 years and then fit an expo¬ 
nential decay function to only the flux measurements of the low- 
energy lines. This yields a 44 Ti mass of (1.28 + 0.14) 10 ~ 4 M G . 
The same fitting procedure applied to the 1157 keV gamma- 
ray line flux measurements (COMPTEL and SPI) results in a 
44 Ti mass of (2.72 ± 0.43) 10 ~ 4 M Q . Using the first value as an 
estimate for the real observable 44 Ti mass seen to decay, that 
is, uncontaminated by possible secondary processes that might 
mimic additional 44 Ti, we evaluated an additional flux level of 
(2.29 ± 0.62) 10 ~ 5 ph cm ~ 2 s _1 by fitting the high-energy mea¬ 
surements with a constrained 44 Ti mass. This discrepancy has a 
statistical significance of 3 cr. 

About 340 years after the explosion, the expected flux ratio 

/F\ 157 is — 1, the measured flux ratio is ~ 0.5. This would 
imply an efficiency (number of photons per decay) error of 50% 
and is excluded by several measurements ([Norman et al. 1998 
IWietfeldt et al.|1999HAhmad et al.|2006) . 


3.3. Nuclear de-excitation? 

We speculate that this discrepancy may be due to an additional 
nuclear de-excitation component originating from accelerated 
particles colliding with 44 Ca* in the remnant and swept up in¬ 
terstellar medium. CRs are thought to be accelerated in SNRs by 
diffusive shoc k acceleration ( Baade & Zwicky 1934), in partic¬ 
ular in Cas A ( |Berezhko & Volk|2004| JAharonian et al.|2001) . 
Only the 1157 keV line from excited 44 Ca will have an addi- 
tional contribution from nuclear excitation (see Sect. \2.\) . The 
additional gamma-ray line flux originating from the isotope 44 Ca 
can be estimated by 


F y = (And 2 y l ■ n( 44 Ca) 


/ 


Qp(p)cr(p)v(p)dp, 


(7) 


where d is the distance to Cas A, n( 44 Ca) the mean density of 
44 Ca atoms in the medium where the interaction takes place, 
Qp(p ) the proton acceleration spectra 110 cr(p) the cross section 
of the reaction 44 Ca(/i, //) 44 Ca*, and v(p) the velocity of accel¬ 
erated protons. The product of n( 44 Ca) by the integral, /( 44 Ca) in 
Eq. ([7]t describes the collision frequency and includes the present 
44 Ca abundance in the interaction region (Ra maty et al.|[l979 
|Summa et al.|201 T) . 

With Eq. ([ 7 ]), the 44 Ca density required in front of the shock 
to produce this gamma-ray flux can be estimated. The mean 
cross sections of the reactions 12 C(p, p') 12 C* by R amaty et al. 
(1979), and 44 Ca(p,//) 44 Ca 


by 


Mitchell et al. 


(1982) have 


the same order of magnitude. With the solar abundance ratio 


1311 


2 ■ 1(T 4 by 


44 Ca/ 12 C = 

gral /( 44 Ca) can be estimated to be 2 


7.001to 6 


Lodders 


10 4 /( 


(2003), the inte- 
12 C). /( 12 C) was 


2 The proton acceleration spectrum, Qp(p) oc p~ 23 , was chosen fol¬ 

lowing Su mma et al.H20lT) to meet the FERMI-LAT constraints on the 

energy content of the cosmic rays that were measured by Abdo et al.| 

(2010), about 3-4 10 49 erg. 


calculated from the values given in Summa et al.]( |2011] >. The re¬ 
sulting 44 Ca density is about 10 5 - 10° cm~ 3 . This value is very 
high compared to densities of the interstellar medium, but may 
reflect the swept-up density in front of the forward shock (“snow 
plough”). If the measurements of presolar grains by Nittler et al. 


( 1996]), Clayton et al. ( 1997| ), and Hoppe et al. ( 20C)0| ) are taken 
into account, which suggest that the ratio 44 Ca/ 4( Ca can be more 
than 100 times higher than solar, the forward-shock 44 Ca density 
is scaled down correspondingly to 10 3 - 10 4 cm -3 . Although this 
density range appears plausible for optical knots and might even 
be consistent with expectations (Fesen et al.||2006[ Docenko & 


Sunyaev [2010j >, the estimated values appear extreme: Assuming 
a volume of the swept-up material near the Cas A SNR of 
~ 10 56 cm 3 , and a density of ~ 10 3 cm~ 3 , the estimated 44 Ca 
mass would be of the order 10 3 M G . Therefore, the ejecta from 
Cas A alone cannot explain this mass; adding the swept-up in¬ 
terstellar material cannot explain the additional flux either. The 
morphology of the swept-up material possibly differs from the 
expectations, or the physics of nuclear excitation may not be 
properly treated in our simple estimate. It would be interesting 
to investigate the nuclear de-excitation lines from 12 C and 16 0 in 
Cas A because these lines are expected to be more luminous. 

Additional support for our speculated nuclear de-excitation 
component is the measured Doppler velocity of the 1157 keV 
line: It is lower than the one derived from the 78 keV line. While 
the emission originating in the decay of 44 Ti to 44 Sc (78 keV 
line) follows the kinematics of the expanding SNR, the 1157 keV 
line emission might have two components, one resembling the 
kinematics (-5000 km s” 1 ) from the 44 Sc decay to 44 Ca, and 
another one incorporating the nuclear de-excitation of 44 Ca* in 
the swept-up material in front of the shock at zero velocity. 
However, we do not know where the interaction takes place and 
can also assume that a measurement of the 1157 keV line sam¬ 
ples another volume element of Cas A than the 78 keV line. 


Milisavljevic & Fesen (2015 (illustrated the bubble-like behavior 


of the interior of Cas A that naturally shows velocities of about 
2000 km s 1 near the expansion center and up to 5000 km s 1 
farther outside. 


4. Conclusions and discussion 

Revisiting the 340 year old supernova remnant Cassiopeia 
A with eleven years of data from the spectrometer SPI on 
INTEGRAL, we detected two gamma-ray lines originating in 
the 44 Ti decay chain 44 Ti—> 44 Sc—> 44 Ca. One of the low-energy 
lines from the process 44 Ti—> 44 Sc at 78 keV was detected, while 
the adjacent line at 68 keV coincides with a major and difficult 
instrumental-background feature and cannot be measured. The 
78 keV line parameters are consistent with other measurements 
constraining the 44 Ti mass and kinematics in the supernova rem¬ 
nant, ((1.5 ±0.4) 10- 4 Mo, and u M ^(4300 ± 1600) km s^ 1 ). The 
second decay step in the 44 Ti decay chain, 44 Sc—> 44 Ca, occurs 
only a few hours after the 44 Ti decay and produces character¬ 
istic 1157 keV gamma-rays that were also detected and ana¬ 
lyzed independently. The 44 Ti mass from the 1157 keV line is 
(2.4 ± 0.9) 10 4 M G , while Doppler broadening here implies an 
expansion velocity of (2200 ± 400) km s~*. 

The kinematic information of the Cas A SNR by mea¬ 
suring emission lines is not unique. In particular, fast-moving 
O-rich knots at velocities of around 8000 km s 1 have been 
found (Fesen et al.] |2006| i. From Chandra measurements, char¬ 
acteristic X-ray lines of Si and Fe have been mapped. An inver¬ 
sion of the classical onion-shell structure had been inferred, as 
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Fe X-ray emission clumps have been found outside the Si X- 
ray emitting regions more central to the remnant (Hwang et al. 


2004 ). From th e NuStar image and the more central location of 
44 Ti clumps (Grefenstette et al. 2014 1 , this inversion hypothe¬ 
sis may be recognized as a misinterpretation of the Fe line X- 
ray emission, which is driven by the ionization state of rem¬ 
nant gas; Fe in the central region thus may escape X-ray line 
detection. The forward shock reached a diameter of 5 arcmin 
and appears to currently move at 5000 km s -1 , while the X-ray 
bright regions predominantly are beyond the inner 2 arcmin, and 
in the jet regions at radial distances above 3 arcmin from the 
central compact object (Hwang et alTj 2004). Thus the reverse 
shock has progressed about half into the remnant (inner diame¬ 
ter ~2-3 arcmin) and has not yet ionized the inner regions (see 
also (Milisavljevic & Fesen 2015 i and (Grefenstette et a l.|2015 i 
for discussions in NIR, and X-ray continuum wavelengths, re¬ 
spectively). 

44 Ti nucleosynthesis yields in core-collapse supernovae is 
one of the questions we wish to address; t hey could range 


up to several 10 5 M Q (Timmes et al. 


1996 


The et al. 


2006 


Magkotsios et al.| 2010|), while for thermonuclear supernovae 


(SN la) only a rare subtype is considered a plausible 44 Ti source, 
but then also with yields of up to several 10 5 M Q (Woosley & 


Weaver 1994), see also The et al. (2006). Nuclear burning in 


supernova explosions occurs at high densities and temperatures, 
and thus those sites are closest to conditions of equilibrium burn¬ 
ing, where nucleons find their most stable arrangement within 
atomic nuclei. The binding energy of nucleons is maximized for 
56 Ni, and nuclear statistical equilibrium (NSE) appears to be a 
good description of plasma conditions for the brief moment of 
the supernova explosion when nucleosynthesis occurs. Here all 
nuclear reactions except weak-force reactions are in thermal bal¬ 
ance. As the explosion site is diluted and cools, the first nuclear 
reactions falling out of such equilibrium are three-body reactions 
such as the 3 a conversion of helium to carbon. Here, the condi¬ 
tions are better described by an overabundance of a particles 
over NSE. This is called a -rich freeze-out, and nucleosynthesis 
productions emerging from these conditions are characterized by 
overabundances of nuclei that are multiples of a particles result¬ 


ing from successive captur es of as on nuclei (Bodansky et al. 
1968 Woosley et al. 1973 1 . 40 Ca and 44 Ti are among the most- 


massive of such products; note that 56 Ni is an a-multiple nucleus 
as well. It has been considered plausible, therefore, that super¬ 
novae that produce major amounts of 56 Ni are also sources of 
significant amounts of 44 Ti. 


When the measurements taken by various instruments dur¬ 
ing the last 20 years are combined, the 44 Ti mass seen to decay 
is (1.37 + 0.19) 1CF 4 M 0 . But we find a significant difference 
(3 cr) in flux, resulting 44 Ti mass, and velocity from the mea¬ 
surements of the 44 Ca versus 44 Sc lines (high- and low-energy 
lines, respectively). We argue that the low-energy line measure¬ 
ments (68 keV and/or 78 keV) reflect the true 44 Ti content and 
kinematics, and an additional component with a flux level of 
(2.29+0.62) 1CF 5 ph cirT 2 s -1 is found in the 1157 keV line flux. 
We speculate that this excess is due to the nuclear de-excitation 
of 44 Ca* in the remnant and swept-up interstellar medium, ex¬ 
cited by cosmic rays that are expected to be accelerated in the 
young remnant by diffusive shock acceleration. 
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